Parameter free calculation of hadronic masses from instantons. 
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We propose a non-perturbative calculation scheme which is based on the semi-classical approx- 
imation of QCD and can be used to evaluate quantities of interest in hadronic physics. As a first 
application, we evaluate the mass of the pion and of the nucleon. Such masses are related to a 
particular combination of Green functions which, in some limit, is dominated by the contribution 
of very small-sized instantons. The size distribution of these pseudo-particles is determined by the 
't Hooft tunneling amplitude formula and therefore our calculation is free from any model param- 
eters. We prove that instanton forces generate a light pion and a nucleon with realistic mass ( 
M„ ~ 970 MeV ). In connection with sum-rules approaches, we discuss the overlap of instantons 
with pion and nucleon resonances. 
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' After the discovery of the BPST instanton solution in QCD [Q, two important results by 't Hooft originated a 
large interest in the semi-classical description of the QCD vacuum ^: the calculation of the tunneling amplitude 
and the discovery of the existence of quark zero-modes in the field of an instanton (which allowed to solve the U(l) 
problem) 0. 

Unfortunately, it was immediately realized that in QCD (unlike in the electroweak case) the tunneling amplitude 
I at one or two-loops accuracy is divergent in the infrared and therefore that instantons "tend to swell" . This fact 
psj ■ led some to the conclusion that the QCD vacuum should be populated by large sized, topologically non trivial field 
CO ' configurations, which cannot be described semi-classically. 

On the other hand, in Q Shuryak argued that the divergence of the single instanton density is non-physical and 
that the size distribution of instantons should eventually start to decay above some average size p, either due to 
non-perturbative effects (which are not included in the 't Hooft formula) or due to interaction with other instantons 
1^ (for a detailed discussion of this issues, see Q). Along this line, he concluded that a semi-classical approach can 
still justified if a small fraction of the Euclidean space time is occupied by instantons, i.e. if p"* n << 1, where n is the 
average instanton density. From an analysis of the global properties of the vacuum (quark and gluon condensate) he 
estimated these parameters to be n ~ 1 /m~^, p ~ 1/3 fm. Based on such arguments, the Instanton Liquid Model 
D (ILM) was developed in the 80 's and 90 's. The model was proven to be phenomenologically successful in describing 
many important low energy properties of QCD, such as the breakingof chiral symmetry Q, several hadronic masses 
^ I and correlators 1^ and, recently, electro-magnetic formfactors [|[ llT| , [l2| . Moreover , studies of the classical content 
of the QCD vacuum, based on lattice simulations via cooling procedure |l6t , confirmed the presence of topologically 
non trivial "bumps" of gauge fields, with typical density and size very close to that assumed by the ILM. 



. Despite its phenomenological success, the ILM remains a model, because none has been able to derive from QCD 
■ " ' why large instantons are not present in the vacuum. However, the fact that there are small-sized instantons in the 
vacuum, whose distribution is described (at two-loop level) by the pAgco << 1 limit of the 't Hooft tunneling 
amplitude, is a result of QCD. 

In this work, we try to set up a calculation scheme for hadronic observables which emphasizes the contribution 
from such small-sized instantons, and we apply it to evaluate the mass of the pion and of the nucleon. In practice, 
in order to keep our calculation finite at each intermediate stage, we will need to regularize the size distribution by 
cutting-off large instantons. However we shall show that, in an appropriate limit to be defined below, at the end of 
the calculation our results will not depend on the particular choice of the cut-off and the only parameter left is Aqcd ■ 

This way, we derive the existence of a light (/ = 1, Jp = 0~) meson state (pion) and of a (/ = i, = i^) baryon 
state (nucleon) . In the case of the nucleon, we are able to extract its mass and find (3.9 ± 0.2) Apv- However, this 
framework did not allow us to determine with sufficient accuracy the very small mass of the pion. 

The paper is organized as follows. In section II we introduce our particular co mbi nation of two-point functions and 



show how this is related to the masses of the corresponding hadrons. In section HI, such a combination is evaluated 
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theoretically, in our semi-classical approach, 
section 



Results are presented and discussed in section [V and summarized in 



II. DETERMINATION OF THE HADRONIC MASSES FROM POINT-TO- WALL CORRELATORS AT 

LARGE MOMENTUM 



Point-to-wall hadronic correlators are defined as the spatial Fourier transform of the point-to-point Euclidean 
correlation fimctions. 



Gf(i,p) = J d^xe'P-^ < 0| JH(t,x)4(0)|0>, 



(1) 



where Jh is an interpolating operator, which excites states with the quantum numbers of the hadron H. These 
correlators are often used to extract hadronic masses from lattice simulations. In fact, by computing their spectral 
decomposition, one finds that, in the large t limit, the contribution from all excitations is exponentially suppressed. 
The energy of the ground state can than be determined from the logarithmic derivative 
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(2) 



By setting p = in one can directly extract the mass of the lowest lying particle, with the quantum number of 
H. In lattice simulations, p is typically chosen to be the smallest possible momentum compatible with the box size 
and dt to be the lattice spacing in the time direction. 

In this work, however, we will need to consider ratio (||), for |p| > 1 GeV. In fact, we shall show that only in such 
kinematic range the pion and nucleon correlators are dominated by the effects of single small-sized instantons, which 
are calculable theoretically. 

As we increase the momentum, the determination of the masses gets more and more difficult, and eventually 
becomes impossible if |p| >> mif. This is due to the fact that, in the ultra-relativistic regime where -E(|p|) ^ |p|, 
all states in the spectral decomposition have the same exponential decay As a consequence, as we increase the 
momentum, larger and larger Euclidean time intervals are needed, in order to isolate the ground-state. Therefore, 
for a given value of the time interval t, the LHS of eq.(^ will be dominated by the ground-state only for a limited 
range of momenta. Generally speaking, we expect the LHS of eq.(^) to become bigger than the RHS, at some large 
momentum. 

The importance of the contribution form resonances depends on the details of the spectral function. In the case 
of the pion, we expect that they should become relevant only at very large momenta, because the first resonance 
(7r(1330)) is much heavier than the ground-state. To estimate such contribution, let us consider a very simple two- 
poles model, in which only tt and 7r(1330) states are included. Than the spectral decomposition of G'2(p,t) would 
read: 



G2{p,t) 



2^? 



OTr(1330) 



0) *) 



2./p2 -hm^ 



(3) 



(1330) 



where and a. 



ir(i330) a-re the coupling of the states with the interpolating operator (which we shall assume to be 
of the same order of magnitude). For t — 7GeV^^ and |p| ~ 2 GeV we have that the contribution of 7r(1330) is 
suppressed by one order of magnitude. 

The nucleon channel is more problematic because the first resonance is only about 500MeV heavier than the 
ground state. In this case, by repeating the above analysis, we find that the desired suppression should be achieved 
for t ^ 8GeV^^ and |p| ~ 1.5GeV Above such thresholds, a theoretical analysis of (||) can be used to obtain 
qualitative information about the strength of the resonances and the continuum. 



III. THEORETICAL CALCULATION 



In this section, we evaluate the LHS of (^) in a semi-classical approximation. Our starting point is the tunneling 
amplitude for one instanton, which is derived by integrating over small quantum fluctuations about the instanton 
solution. At two-loop accuracy it reads: 



U-'tHooft 



const ■ p ^ p^' 



Pi{p) 



2N^ 



exp 



-P2{p) 



b' 

2Nc 

2b 



2&/3i(p) 



ln[/3,{p)] 



(4) 



3 



where f3i{p) and /32(p) are the one and two- loop beta functions: 

Piip) = -b In(pApv), 



/32(p)=/3i(p) + |^/n(^/3i(p)), 



(5) 
(6) 
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(7) 



This resuh was carried out in the Pauh-Villars regularization scheme and Apy is the corresponding scale parameter. 

This formula is derived perturbatively and therefore is valid only for small instantons, for which pAgcD << 1- 
Its infrared divergence indicates that a fully perturbative description of a single tunneling amplitude is inconsistent. 
This is the so called 'infrared catastrophe' of the instanton calculus in QCD. 

In the ILM, this problem is overcome by assuming that the vacuum can be described as an ensemble of interact- 
ing instantons, and that the single instanton size distribution is cut from above by an average instanton-instanton 
repulsion. The model depends on one dimensionless parameter, the instanton packing fraction p'^ n, which can not be 
derived and has to be estimated phenomenologically. 

Such an approach necessarily introduces some additional approximations. First of all, in the presence of non 
negligible overlapping, the sum of many instantons is no longer an exact solution of the classical equation of motion. 
Moreover, the tunneling events can non be treated as completely independent and therefore it is not guaranteed a 
priori that the contributions of such classical configurations is not overwhelmed by quantum fluctuations. However, 
one should mention that the ILM gives indeed a very good phenomenology, therefore one can argue, a posteriori^ that 
the assumptions made are legitimate. 

Here, we try to go one step further and set up a calculation scheme which is sensitive only to the small-size tail of 
the instanton density and contains no phenomenological parameters. We shall only assume that the instanton packing 
fraction is small enough for the pseudo-particles not to loose completely their individuality. In other words, the QCD 
vacuum will still be described as some instanton ensemble. However, we will not need to choose a particular value 
for p and n, nor we will have to model the instanton-instanton interaction. As a matter of fact, we shall also not 
require that the semi-classical approach is justified for all pseudo-particle in the ensemble, because our results will be 
sensitive only to the contribution from the small ones. 

The basic simplification introduced by the semi-classical approximation is that the infinite dimensional integration 
over all gauge configurations is replaced by a finite number of integrals over the collective coordinates of each pseudo- 
particle. 

Let us consider, for sake of definiteness Q, the pseudo-scalar two-point correlation function, interpolated by the 
current: 



J5(x) = u{x)j5d{x). 

In our approach, the two-point correlator is approximated with (for two flavors): 



(8) 



G{x,y) :=<0| J5(a;)4(0)|0>~| J dnf{n) 



Yl det{IJ)n + mf) 



Tr[-fr,S[x,Q-n)-f^S{Q,x;n% 



(9) 



z = j dnfivt) 



W det{I?n+mf) 



f= 



(10) 



where = (fii, f2jv+ , ^^Ar++i, ■■■,^Ni+Na) is the set of all collective coordinate identifying the configurations of the 
instanton ensemble /(O) is the weight associated to a particular configuration in the pure Yang-Mills theory 
p3| and S{x, y; fl) is the quark propagator in the configuration il. The quark propagator can be formally written as: 



S{x,y;n)^Y. 



—i X + m 



(11) 
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where V'a (■^) eigenvalues of the Dirac operator, TjP'^^ = ^^'^ field of a single instanton, the propagator 

is known exactly [[l4| and consists of a non zero-mode part and a zero-mode part: 



S{x,y; Qi) = Snzm{x,y; Qi) + Szm.{x, y\ fii), 



(12) 



where J7i denotes the set of collective coordinate relative to the particular instanton i. It can be verified that, when 
correlators receive contribution from more that one zero-mode propagator ( like in the case we are considering), the 
non zero-mode part can be very well approximated with the free propagator: 



(13) 



The propagator in the field of many instantons was evaluated in and consists of a non zero-mode part (which can 
again be approximated with a free propagator) and a quasi zero- mode part, which can be expanded in the functional 
space of zero- modes of individual instantons [p4[ : 



(14) 



I, J 



where ip^' (x) are the well known zero-mode wave functions ||^, | 
Inserting ( p^ in (H) we find: 



N++N- N++N- 



G2{x,y) 



G2free{x,y) = ^ I J2 E / ^^/(") 

V / j=i /' j'=i 



Yl det{I])n + mf) 

f=u,d 



r J' 



■Tr[j, ^P^' (x) V?'^ \0) 75 ^o'^' (0) 4'o-" \^)] 



(15) 



where N^^ ^ denotes the total number of instantons (anti-instantons). It is convenient to re- write this sum in a such 
a way that the contribution of a single instanton (anti-instanton) is isolated: 



G2 (X,y)~ G2 free {x , v) = ^ \ N+ j d nf{n) 



Yl det{I^n + mf) 



C?,fTr[j, Vo ^ {x) ^ ^(0) 75 V'"^ (0) i''o''\x)] + 



-N- J dflfifl) 



Y[ det{pn + mf] 

f=u,d 



IC-.+i A..+irrr[75 Vo'"""^ ix) ^^"'''\0) 75 V'o (0) V^""' ' (^)] 



^N + + i /n^ + + i t / 



(16) 



where the ellipses denote terms which depend on the collective coordinates of more than one (anti)instanton. Such 
terms will introduce effects which are suppressed by powers of the instanton density and are typically sub- leading, as 
will be discussed later. Now we can integrate out the remaining + — 1 components of and find: 



1 



G2ix,y)-G2free{x,y)^{ / d f^i /(f^i ) — Tr [75 V'o ^ (^) V^o ' '(0) 75 V'^ (») ^0 H^^)] + 
+ [ rff^^++i/(f^Ar++i) ^Tr[75 V?"^"' W V'o "^^^ ^(0) 75 V'o "^^^ (0) i^o^^^'hx)] 



where we have introduced the single-instanton distribution p£ 



(17) 



f{ni)^N+/z dn2...dnN+f{n) 



Yl det{I^n+mf) 



f=u,d 



(18) 



and the mass parameter m* has been added to keep track of the dimensionality. 

Now let us now discuss the further simplifications that arise by applying just symmetry principles. First of all, 
since the correlator we are considering is gauge invariant, the integration over the color orientation of the instanton 
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very long point-to-point correlator 
for whiclr the SIA does not apply. 




point-to-point correlator 
for which the SIA applies 



FIG. 1: Contribution of different point-to-point correlators. For p ~ IGeV only tlie correlators with one end-point in the 
shaded area contribute. For such correlators, the SIA is supposed to hold. 



and of the anti-instanton is completely trivial. Moreover, total momentum conservation requires that the distribution 
of such instanton in the vacuum is homogeneous |2q|, that is: 



f{z;p) ^nd{p) 

Notice that we have completely recovered the structure the single instanton density introduced by Shuryak in 
proposed the very simple ansatz: 

f{z;p) ^n6{p-p), 



(19) 
who 

(20) 



which emphasizes the global properties of the instanton vacuum (i.e. average instanton size and density). 

In our calculation, however, we want to focus on the effect of small instantons. At this purpose, we choose a 
parameterization that, in the limit p — s- 0, reproduces the tunneling amplitude (^. For example we can choose the 
following regularization: 



1 



1 + exp{{p - p)/a) 



(21) 



We can now rewrite (|T^ as: 



Gl'^ix^y) = G2{x,y) -G2free{x,y) = n[ j d*z j d pd{p)Tr[-fz Sl^{x,0; z, p)-/^ Sl^{0; x; z, p)] + 

d^z I dpd{p)Tr[-/5S^^{x,0;z,p)-f5S^^{0;x;z,p)] 



(22) 



where we have taken equal density of instantons and anti- instantons, nj = ua ='■ n, and Szm {x, y; z) is the zero-mode 
propagator of a particle of mass m* in the field of an instanton (anti-instanton) of size p at z. Hence, we can identify 
the parameter m* as the quark effective mass originally introduced in ||^ . In |^] such parameter was estimated from 
numerical simulations and form phenomenology to be m* ~ S'oMeV . 

Summarizing, what we have done so far is to provide a formal derivation of the Single Instanton Approximation 
(SIA), starting from the partition function (px|). The validity of such approximation, in the case of the ILM, was 
studied in detail in |p^ . It that analysis, we considered several point-to-point correlators that receive contribution 
from two zero- mode propagators and compared SIA calculations (with the ansatz (p9|)) with numerical simulations 
in the Random and Interacting Instanton Liquid. We found complete agreement up to distances of the order of 
6 — 7 GeV~^ , depending on the particular correlator. Since this value was obtained in a particular model, it should be 
considered as a rough estimate. We may argue that, if a particular observable is mainly sensitive to the contribution of 
small instantons, the SIA should somehow increase its range of validity, because the packing fraction of such pseudo- 
particles is smaller. Also, we do not expect such SIA calculation for the point-to-wall propagator to be reliable for 
all values of the momentum p. In fact, if such momentum is small, the Fourier transform (y) receives non-negligible 
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t/2 t time 

FIG. 2: Interpretation of the single small instanton domination of the point-to-wall correlators. 

contribution from point-to-point correlators which connect the origin with points on the wall that are very far away 
from the time axis (see fig. |^). From the discussion above it follows that such point-to-point correlators have size 
larger than the maximal compatible with the SIA. However, for |p| of the order of a GeV, only points at the distance 
smaller than roughly one inverse GeV from the time axis will contribute to the Fourier transform, and the SIA is 
applicable. 

Summarizing, the feasibility of our approach relies on the existence of a window in the time and momentum 
variables, such that the SIA is still holding and the ground-state is not overwhelmed by resonances. For example, for 
t > 9 GeV^-^ and p < IGeV, we do not trust our approximation, while for t < 6 GeV^^ and p ^ IGeV, the SIA would 
be reliable, but the correlation functions would start to become very sensitive to the contribution of the resonances. 
Based on such arguments, we shall choose t to be 7 GeV^ ^ t ^ SGeV"^ and |p| ^ 1 GeV. 

Once the general calculation scheme has been established, we can go ahead and evaluate the LHS of (||) for pion and 
nucleon. We begin with the pion pseudo-scalar point-to-wall correlation function. The leading instanton contribution 
to (i, p) is given by: 



where f (x) = + . 

The stationary phase method can also be applied to determine what instantons contribute the most, as |p| increases. 
By expanding the exponent around its minimum (p = 0, Z4 = t/2), it is straightforward to verify that, at large 
momenta, G2(t, |p|) will be dominated by small-sized instantons. Hence we expect that, in this regime, the prediction 
should be fairly insensitive to the detail of the cut-off function in and depend essentially on nitHooft, only. 

The result that a single, small instanton dominates the point-to-wall correlation function at large momentum has 
an intuitive interpretation p7| . The current at origin emits a quark and an anti-quark with undetermined momenta 
(see fig. H). In order for such particles to form a bound-state at some later time t, their momenta has to be focused 
and made almost collinear. The nearby instanton works as a "lens" and provides the required bending of the "quark 
beams" . At high momentum, the focusing is optimal if the "curvature" of the instanton lens is small and sits at a 
distance 2:4 = 1/2 1 from the source. 

The same calculation can be repeated for the nucleon point-to-wall correlation function. In such case, we have 
more freedom in the choice of the point-to-point correlator. Previous studies Q have shown that the single-instanton 
contribution is largest in the correlator: 

G'iix, v) Q\Tr[3Ux)~jsM n\0 >, (24) 

where jsc{x) = [ua{x) {G df,(x)\uc{x) eabc is the so called 'scalar current'. The leading instanton contribution to 
the nucleon point-to-wall correlator reads: 

Gr ft,P) = 11 ,22 / ^ ' ; / dpd{p)p^ dzi d|z| 



pK2^2y^ p+i2]3/2y^ 7o ' '[|zP+e(^4)]2[(r2 + |z|2+e(z4-t))2 



(25) 
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Again, we expect this result to be reliable only for |p| ^ 1 GeV and that small instantons should dominate, in the 
large momentum regime. 

One can check that, at the scale we are working at, the free contribution to the pion and nucleon point-to-wall 
correlators is much smaller than the corresponding leading instanton term, for any reasonable value of m* . Hence, 
in the ratio (|^), the free contribution can be safely neglected and the model parameters n and m* cancel out. As 
a result, the only parameters that remain are Aqcd, appearing in d'tHooft{p) and the position of the instanton size 
cut-off p. However, from the above discussion, we expect the results to become independent on p for sufficiently large 
values of p. 

In this section we have focused on the evaluation of the logarithmic derivative (||) , which allows to extract hadronic 
masses. In principle, the same calculation scheme can be applied to other hadronic observables such as, e.g., electro- 
magnetic form factors . Following this prescription, one needs to relate the observable being investigated to some 
ratio of point-to-wall correlation functions and choose t and p according to the above discussion. One should keep in 
mind, however, that the SIA approximation is much more accurate if the relevant Green functions receive contribution 
from more than one quark zero-mode propagator |17|. 



IV. RESULTS AND DISCUSSION 



p=1.6 GcV 
p = 2.0 GeV" 
p = 2.5 GeV "' 




P [GeV'l 

FIG. 3: The (unnormalized) single instanton tunneling amplitude for different values of the cut-off p {A.pv = 250 MeV). 



Let us begin by discussing the results obtained in the case of the pion . In fig. ^, several theoretical evaluations 
of the pion dispersion curve i?(|p|), are compared with the physical one. The different SIA predictions were obtained 
from (g) using different values of the time interval and of the size cut-off |^ see fig. ||). Notice that these values 
correspond to extremely different pictures of the vacuum, because the corresponding packing fraction is proportional 
to p^ and changes by one order of magnitude. 

First of all we can verify that the pion pole has been isolated and that the results are fairly independent on the 
value of cut-off used. From the comparison in the region IGeV ^ |p| ^ 2GeV, we conclude that our semi-classical 
approach provides with the right amount of dynamics and is completely consistent with the existence of a light pion. 

Following the discussion of section II, we would expect the excited states to contribute to the relevant green- 
functions, in the region above 2 GeV. This implies that the logarithmic derivative (H) must deviate from the pion 
dispersion curve, at some high momentum. However, we observe that the complete agreement between our calculation 
and the pion dispersion curve continues up to high momenta. This fact can be interpreted in terms of quark-hadron 
duality p8| , which assumes that the contribution from the continuum is dominated by perturbative effects, while 
the ground-state pole can be described in terms of non-perturbative effects. Indeed our purely non-perturbative 
calculation does not reproduce the presence of a continuum. 

The results obtained in the case of the nucleon, using different Euclidean time intervals and instanton size cut-offs, 
are reported in fig. ^. The numerical solution of the integral (25) is quite challenging and the signal becomes very 
noisy above 2 GeV. Again, we verify that the ground state pole has been isolated and that the results are insensitive 
to the choice of the cut-off. We observe agreement between our theoretical calculation and the experimental curve in 
the expected region, IGeV < |p| ^ l.GGeV^ . Moreover, unlike in the case of the pion, the mass of the nucleon is 
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FIG. 5: Dispersion curve of the nucleon {Apv = 250 Me V). The sohd hne is obtained from the physical nucleon mass, the 
points are SIA results. 



sufficiently large do be extracted from the best fit curve |30 . We found M„ = 3.9 ± 0.2Apv' ~ 970MeV, consistent 
with the experimental value. 

In this case, by extrapolating the trend of the calculations above 2 GeV, we may say that some little deviation from 
the nucleon dispersion curve seems to occur, at high momenta. This fact may suggest that, unlike in the case of the 
pion, instanton non-perturbative effects have some overlap with nucleon resonances. 



V. CONCLUSIONS AND OUTLOOK 



In this paper we have presented a new calculation scheme which can be used to study several hadronic observables. 
The method is based on the idea that, in the appropriate limit, it is possible to isolate the contribution form small-sized 
instantons, upon which we have theoretical control. 

The prescription was applied to evaluate the mass of the pion and of the proton. We found agreement with 
phenomenology and verify that our results are insensitive to the particular choice of cut-off p and therefore are model 
independent. 

By comparing the results obtained for the pion and nucleon, we discussed the importance of instanton effects for the 
physics of resonances, in these channels. We found that, in the case of the pion, a purely non-perturbative instanton 
calculation does not reproduce the existence any resonance or of a continuum up to extremely high momenta. In the 
case of the nucleon, some deviation from the ground-state dispersion curve seems to occur at high momenta (although 
our results are not accurate enough to make a definitive statement). This would imply that instanton effects could 
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have some overlap with nucleon resonances. We intend to develop a more quantitative analysis of this issue in a future 
work. 

We think that the method could be useful in investigating the contribution of instantons to a variety of hadronic 
observables. An application of the same ideas to the analysis of pion and proton electro-magnetic formfactors is in 
preparation. 
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The definition of the single anti-instanton is, of course, completely equivalent. 

The implication of total momentum conservation in SIA will be discussed in detail in ]l^ . 

We would like to thank E.V. Shuryak for suggesting such interpretation. 

In performing the calculation, we have expressed all dimensional quantities in terms of Apv-The results presented in this 
section have been converted to natural units by choo sing Kpv ~ 250MeV. 

These results were obtained using a = O.lGeV^"^ in (pl|). Different choices were tried and gave equivalent results. 
According to the discussion in (II), we have fitted only the points in the region 1 GeV < |p| < 1.6 GeV. 



